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THE  PROBLEM 


In  the  formulation  of  lumped  parameter  models  to  predict  the  dynamic  response 
of  the  human  body  such  rigid  body  characteristics  as  constant  mass,  center  of  mass, 
and  moments  of  Inertia  are  ascribed  to  anatomical  segments.  The  subject  of  this 
report  Is  a technique  for  direct  measurement  of  these  quantities. 

FINDINGS 

The  measurement  techniques  were  performed  on  embalmed  specimen  human 
heads  and  human  head  and  necks.  The  mass  distribution  parameters  as  measured 
are  within  five  percent  of  the  actual  specimen  value,  but  the  specimens  themselves 
are  changed  considerably  from  the  living  state  as  a result  of  the  embalming  procedures. 

RECOMMENDATIONS 

These  measurement  techniques  can  be  successfully  applied  to  anatomical  segments 
of  Interest,  but  care  should  be  taken  to  keep  the  specimen  segment  from  deteriorating. 
The  most  promising  method  of  preservation  at  this  time  Is  the  deep  freezing  of  the 
anatomical  segment  as  soon  as  possible  In  order  to  stabilize  the  body  fluids  within  the 
segment. 
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Measurement  of  Mass  Distribution 
Parameters  of  Anatomical  Segments 

Edward  B.  Becker 

Naval  Aerospace  Medical  Research  Laboratory 
Abstract 

Procedures  to  determine  the  center  of  mass  and  the  moments  of  inertia  in  three 
dimensions  of  previously  defined  anatomical  segments  are  presented.  As  an  illus- 
tration, these  procedures  are  applied  to  the  human  head  and  head-und-neck.  The 
results  of  measurements  made  on  six  human  heads  and  three  hcud-and-necks  are 
presented  and  discussed. 

IN  THE  FORMULATION  of  lumped  parameter  models  to  predict  thedynumic 
response  of  the  human  body,  such  rigid-body  characteristics  as  constant  muss, 
center  of  mass,  and  moments  of  inertia  are  ascribed  to  various  anatomical  seg- 
ments, Since  most  of  these  segments  are  notoriously  not  rigid  bodies,  the  best 
rigid-body  approximation  for  uny  particular  response  could  be  derived  only  from 
an  extensive  study  of  that  response.  Even  so,  the  mass  distribution  of  a particular 
body  segment  for  the  static  case  (that  is,  for  motions  so  slight  us  to  cuusc  no 
deformation  of  the  soft  tissues)  can  serve  as  an  iuitiul  estlmute  of  the  best  rigid- 
body  approximation  of  this  segment  for  mere  complex  behavior. 

The  procedures  outlined  in  this  paper  are  discussed  with  respect  to  a particular 
anatomical  body,  the  human  head.  It  is  felt  that  the  exhaustive  treatment  given 
the  problem  of  upplying  the  subject  measurement  procedures  to  the  humun  head 
und  the  humun  head-und-neck  will  sufficiently  illustrate  the  techniques  for  applica- 
tion to  other  anatomical  sections. 

Any  attempt  to  measure  directly  the  rigid-body  characteristics  of  the  head  and 
of  the  head-und-neck  must  include  a treatment  of  certain  factors,  such  us: 

1.  The  humun  head  und  head-und-neck  have  no  definite  physical  demarcation. 

2.  Neither  the  head  nor  the  heud-and  neck  is  a rigid  body. 

3.  The  physical  structure  of  these  bodies  varies  from  individual  to  individual, 
obscuring  the  selection  of  anatomical  landmarks  with  which  to  loeute  and  align 
these  parameters. 

These  factors  are  treated  us  follows: 

1 . The  physical  limits  of  the  head  and  the  heud-and-neck  are  defined  by  the 
same  unutomicul  procedures  us  those  used  in  Ref.  I . 

2.  Two  factors  compromise  the  rigidity  of  the  these  bodies.  Each  is  composed 
of  jointed  bodies,  which  can  assume  a number  of  positions,  und  each  is  composed 
of  large  amounts  of  soft  tissue.  These  bodies  are  made  to  approximate  rigid  bodies 
by  selecting  a single  configuration  (jaw  shut  and  the  cervical  spine  in  a rigid  but 
uncontrolled  set)  and  employing  measuring  procedures  that  minimize  relative 
motions  of  the  soft  tissues. 

3.  The  auditory  meatuses  und  the  orbital  notches,  which  define  Reade’s  plane, 
are  used  to  define  a three-dimensional  coordinate  system  (Fig.  1 ) in  which  these 
parameters  are  then  loeuted.  In  this  manner,  these  mechanical  characteristics 
are  fixed  in  the  bony  structure  of  the  skull. 
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The  measuring  procedures  employed  require  that  a further  bargain  with  physical 
reality  be  struck  before  any  measurements  are  made.  In  order  to  make  these 
measurements,  the  head  or  the  head-and-neck  must  be  physically  separated  from 
the  rest  of  the  body.  Therefore,  these  procedures  are  performed  upon  cadaveric 
material  in  hopes  that  a statistical  evaluation  of  the  information  obtained  will 
relate  these  mechanical  parameters  to  such  measurements  us  may  be  made  on 
living  subjects. 

Since  the  subjects  on  which  these  measurements  were  made  were  prepared  at 
Tulune  University  for  measurements  of  these  parameters  in  two  dimensions  (that 
is,  the  locution  of  the  center  of  mass  in  the  midsagittul  plane  and  the  moment  of 
inertia  about  un  axis  perpendicular  to  this  plane),  this  study  will  be  concerned 
primarily  with  the  procedures  und  mathematics  ncccssury  to  obtain  the  complete 
center  of  muss  vector  and  moment  of  inertia  tensor.  The  anutomical  preparation, 
the  ages,  und  the  histories  of  the  subjects  arc  reported  in  Ref,  I . 

These  measurement  procedures  are  designed  to  maximize  operational  simplicity 
and  redundancy  ut  the  expense  of  computational  complexity,  In  this  manner, 
procedural  errors  are  minimized,  while  transcription  errors  may  be  recognized 
und  eliminated  during  the  data  reduction  without  compromising  the  results, 

The  bulk  of  these  procedures  involve  measurements  made  on  a single  piece  of 
equipment,  the  stereotaxic  jig.  This  jig  is  composed  of  a David  Kopflndustrics 
Model  1630  stereotaxic  unit  set  into  u tetrahedral  frume.  The  frame  is  designed  to 
facilitate  those  measurements  necessary  to  obtain  the  center  of  muss  und  the 


Fig,  1 • Head-referenced  coordinates 
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moment  of  inertia.  The  stereotaxic  unit  is  designed  to  locate  and  secure  the  subject 
within  it. 

Once  obtained,  the  center  of  mass  and  the  moment  of  inertia  of  the  subject-jig 
combination  and  the  empty  jig  cun  be  compared  to  yield  the  subject  center  of  mass 
and  moment  of  inertia.  Once  the  subject  has  been  located  within  the  jig,  these 
mechanical  parameters  can  be  transformed  to  anatomically  bused  coordinates. 

Procedural  Apparatus 

The  apparatus  cun  be  separated  into  three  groups  by  function:  the  subject  pro- 
cedure interface,  the  center  of  mass  equipment,  it  nd  the  moment  of  inertia  equip- 
ment, 

The  subject  procedure  interface  consists  of  the  jig  and  a DKl  1760  micro- 
manipulator. The  jig  (Fig.  2)  is  a specially  adapted  DKI  Model  1630  stereotaxic 
unit  set  into  a tetrahedral  frame,  This  stcreota.  Sc  unit  secures  the  subject  by  means 
of  two  bars  inserted  into  the  auditory  meatuses,  ,wo  more  bars  plueed  on  the 
inferior  orbital  ridges  over  the  orbital  notches,  and  a fifth  bur  placed  under  the 
chin. 

The  tetruhedrul  frame  is  specially  designed  to  contain  the  stereotaxic  unit  and  to 
facilitate  the  functions  of  the  two  other  equipment  groups.  It  is  built  in  the  form 
of  a tetrahedron  with  projections  at  the  vertexes  and  at  the  midpoints  of  each  edge. 

The  micromunipulator  is  a multijointed  arm  designed  to  position  probes,  elec- 
trodes, and  the  like.  It  is  used  us  outlined  in  Appendix  D to  locate  the  auditory 
meutuses  and  the  notch  on  each  of  the  inferior  orbital  ridges  of  a subject  secured 
in  the  jig.  The  locations  of  these  four  points  then  yield  the  precise  position  and 
alignment  of  the  subject  with  respect  to  the  jig.  The  center  of  mass  equipment 
consists  of  a single-arm  balance  calibrated  to  0.02  lb. 


Fig.  2 - Stereotaxic  unit  and 
tetrahedral  frame 
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By  positioning  the  jig  so  that  it  is  supported  by  three  of  the  edge-midpoint  pro- 
jections and  resting  one  of  the  three  projections  on  the  balance,  it  is  possible  to 
measure  the  load  supported  by  that  projection. 

In  actual  fact,  a problem  intrudes  here.  For  this  particular  equipment  configu- 
ration, there  is  no  stable  equilibrium;  the  baluncc  arm  cun  never  be  leveled,  but 
will  always  be  at  some  extreme  limit,  high  or  low,  of  its  travel.  An  explanation  of 
this  phenomenon  is  included  in  Appendix  A.  The  measurement  made  is  that 
weight  ut  which  the  balance  pan  moves  from  its  highest  to  its  lowest  position,  This 
procedure  introduces  definite  errors  into  the  measurements,  but  these  errors  negate 
each  other  in  the  data  reduction, 

There  are  four  different  positions  in  which  the  jig  may  be  placed  so  that  it  is 
supported  by  three  of  the  edge-midpoint  projections,  yielding  a total  of  12  mea- 
surements that  can  be  made.  These  1 2 measurements  then  yield  the  three- 
dimensional  position  of  the  c.g.  of  the  jig. 

The  moment  of  inertia  equipment  consists  of  three  wires,  whose  ends  may  be 
looped  about  three  of  the  jig’s  projections  suspending  it  in  the  manner  of  a trililar 
pendulum,  and  the  timing  equipment.  The  timing  equipment  is  composed  of  a 
500  W slide  projector,  a mirror,  a photosensitive  trigger,  and  a countertimer.  The 
projector  beam,  trimmed  by  placing  a razor  in  the  slide  receptacle,  reflects  from  a 
mirror  attached  to  the  jig,  and  is  focused  onto  the  photosensitive  trigger.  This 
trigger  is  pluced  us  near  the  projector  us  possible  so  that  while  rotutionul  motion 
of  the  jig  sweeps  the  beam  across  the  trigger,  translational  motion  produces  no 
discernible  beam  displacement.  This  trigger  then  activates  the  countertiiner,  per- 
mitting accurate  measurement  of  the  rotational  period  of  the  jig. 

There  are  four  different  orientations  in  which  the  jig  may  be  suspended  by  three 
vertex  projections,  and  six  more  orientations  in  which  the  jig  may  be  suspended  by 
two  vertex  projections  and  the  edge-midpoint  projection  between  the  two  unused 


Fig.  3 - Moment  of  inertia  procedures 
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vertexes,  as  shown  in  Fig.  3.  The  rotational  periods  of  the  trifilar  pendulums 
formed  by  hanging  the  jig  in  these  10  different  orientations  are  measured,  and 
these  10  measurements  yield  the  full  three-dimensional  inertia  tensor. 


Mathematics 

Most  of  this  section  is  devoted  to  describing  transformations  between  four 
separate  coordinate  systems.  The  first  of  these  coordinate  systems,  shown  in 
Fig,  1,  is  the  one  fixed  in  the  bony  structure  of  the  human  skull. 

Four  anatomical  landmarks,  the  right  and  left  auditory  meatuses  and  the 
notches  on  the  right  and  left  inferior  orbital  ridges,  locate  the  head  coordinates 
using  the  algorithms  in  Appendix  C.  The  X-Y  plane  is  the  plane  that  passes 
through  the  midpoints  of  the  lines  connecting  the  two  auditory  meatuses,  the  two 
orbital  notches,  the  right  auditory  meatuses  and  the  right  orbital  notch,  and  the 
left  auditory  meatus  to  the  left  orbital  notch, 

Proof  thut  these  four  points  arecoplanar  is  provided  in  Appendix  B.  The  first 
of  these  midpoints  is  the  origin.  This  point  and  the  second  midpoint  locate  the 
X-uxis.  The  Y -axis  is  parallel  to  a line  through  the  third  and  the  lust  of  these  mid- 
points, The  X-axis  is  positive  in  the  forward-facing  direction,  the  Y -uxis  is  positive 
on  the  left  side,  and  the  Z-axis,  defined  by  A y,  is  positive  through  the  top  of  the 
skull, 

The  second  of  these  coordinate  systems  is  located  in  the  jig  hurdware.  This  sys- 
tem is  chosen  for  its  rough  alignment  with  the  coordinate  system  of  a subject 
secured  within  it,  and  is  described  in  terms  of  the  subject,  The  origin  is  the  vertex 
behind  and  to  the  right  of  the  right  car,  The  Y -axis  is  along  the  edge  connecting 
the  origin  to  the  vertex  behind  and  to  the  left  of  the  left  ear.  The  X-uxis  is  parallel 
to  the  line  connecting  the  midpoint  of  the  Y-axis  edge  to  the  midpoint  of  the  edge 
on  the  opposite  side  of  the  tetrahedron.  This  edge  is  the  only  one  of  the  five  re- 
maining edges  with  which  the  Y-usis  edge  has  no  direct  contact.  This  particular 
edge  is  also  parallel  to  the  Z-axis  Positive  directions  are  head  left,  up,  and  facing, 
respectively. 

The  remaining  two  coordinate  systems  are  fixed  in  the  procedures.  The  A-  - 
plane  in  the  c.g.  procedures  is  the  plane  through  the  points  on  which  the  jig’s  three 
supporting  projections  rest.  The  origin  is  the  point  of  contact  between  the  pro- 
jection whose  loading  is  being  measured  and  the  balance  pun.  The  A-axis  passes 
through  the  origin  and  is  perpendicular  to  the  line  through  the  other  two  contact 
points.  The  2-axis  is  positive  upward  and  parallel  to  the  gravity  vector. 

The  A - f plane  in  the  moment  of  inertia  procedures  pusses  through  the  loops 
at  the  ends  of  the  three  wires.  The  2-axis  is  positive  upward  and  parallel  to  the 
gravity  vector.  Taking  one  of  these  loops  as  the  origin,  the  P-axis  is  positive 
through  the  next  loop  going  clockwise  (looking  down).  The  A-uxis  is  defined  by 
| x 2. 

In  performing  u complete  set  of  measurements,  the  subject  assumes  one  orienta- 
tion relative  to  the  jig.  The  jig  assumes  1 2 orientations  to  the  c.g.  procedures  and 
10  orientations  to  the  moment  of  inertia  procedures,  each  of  these  22  procedural 
orientations  being  assumed  twice. 

The  orientation  of  the  head  relative  to  the  jig  is  defined  b>  a cosine  matrix  and  a 
vector 
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where: 
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W,  - measured  value  for  ith  orientation 

In  the  moment  of  inertia  procedures,  the  theoretical  value  of  the  pci  iod  of  the 
rotation,  i motion  of  the  trifilar  pendulum  is 


T = 
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as  developed  in  Appendix  F.  or 
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The  value  of /u  in  these  measurements  is  a function  of  the  jig  inertial  tensor 
and  the  orientation. 
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Since /2I  - / ij./ji  - I\}Jd  « /32.  This  equation  cun  be  combined  with  the 
equation  for  the  period  to  yield 

T2„  = <'1(l/||  + cs,,/ 12  + rj„/,i  + c4lt  1 n + c<,„I  ji  + (13) 

where  C|„  - e6„  are  functions  of  the  system  geometry  and  the  location  of  thec.g. 

A least-squares  solution  is  complicated  by  the  fact  that  T„  is  the  quantity 
measured  while  T\  is  the  quantity  mobilized  in  these  equations.  However,  taking 
advantage  of 

(Tj  theoretical  - T2„  measured)2  = (T„  theoretical  - T„  measured)2 

(7’„  theoretical  + T„  measured)2 
Sr  (7'„  theoretical  - T„  measured)2  4 T2„  measured 


and  minimizing  the  term 
yields 

", 


(7’J  theoretical  - Tl„  measured)2/  47', ; measured 

(14) 


(r'l«)  v (*lM*2/t) 


T2 

• it 


J 


n 1 In 


(15) 


.!■  i,  vi  1 1 A HiVf iVJ  .!jrfii&.|S).,u. . 


_ - IvV,  


167 


PARAMETERS  OE  ANATOMICAL  SEGMENTS 


The  results  of  the  least-squares  procedures  for  both  the  e.g.  vector  and  the 
moment  of  inertia  tensor  can  then  be  used  to  calculate  theoretical  values  for  each 
of  the  quantities  measured.  Comparing  these  will  show  which  measurements,  if 
any,  involve  procedural  or  transeriptive  errors, 

These  measurements  are  then  eliminated  from  the  summations  in  a second  ap- 
plication of  the  least-squares  procedures. 

Once  these  procedures  have  yielded  the  inertia  tensor  and  the  e.g.  vector  for 
the  empty  jig  and  for  the  subject  jig  combination,  all  that  remains  is  tocxtruct 
these  quantities  for  the  subject  alone, 


The  subject  e.g.  is  given  by 


I wa)xc.t,.»l  ~ “7-Wl 


similarly 


(16) 


The  subject  inertia  tensor  is  given  by 

f|u  * fi \ij  ~ fit j "E  - ( y c,g.  ij  + eCig>  tj)  — wJ(<yv>g>  j + zUig>  3) 

- wAy^2j  - rv.s.y)2t 

7|2j  “ ^ 1 7sj  ~ I\IJ  ~ l,vjy(-*fc.g.v^’c.g.i/)  — a’, (.Vcg  , t’cg., ) 

- ^(■»'c.s./>’c.g.y)l 


(17) 

(18) 


and  so  on. 

The  values  are  expressed  in  the  jig  coordinate  system  and  ure  transformed  to  the 
subject  coordinates  as  follows: 


“ ’ 

7 

•*<;.g.i 

\ 

»« 

^C.g.i 

- 

*c.g  .s 

subject 

\ 

zc.g.J 

jig 

b3 

/ 

L.  — J 

coordinatr-; 

L J 

\ 

L-  . J 

coordinate* 

L,  _J 

/ 

1 liubject  = !*(()  (Ailjlgl*y) 


where  itJ  and  bt  arc  the  transformation  parameters  in  Eq.  1 . 


(19) 


(20) 


Calibration 

A number  of  small  weights  were  attached  to  the  jig  to  simulate  a body  of  known 
e.g.  and  moment  of  inertia.  These  quantities  were  then  measured  using  the  pro- 
cedures described  in  this  paper. 

The  results  of  the  moment  of  inertia  procedures  were  then  used  to  determine 
empirically  an  operational  value  for  p the  length  of  the  wires  in  the  trililar 
pendulum. 

The  e.g.  of  the  weight-jig  system  and  of  the  empty  jig  was  calculated,  and  the 
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c.g.  was  used  to  calculate  theoretical  values  for  the  weights  measured  at  the  projec- 
tions. As  expected,  these  theoretical  values  differed  in  a uniform  manner  from  the 
measured  quantities.  The  uniformity  of  the  difference  and  the  symmetry  of  the 
experiments  lead  one  to  believe  that  these  errors  are  self-canceling, 

Similarly,  the  moment  of  inertia  tensors  of  these  systems  were  calculated  and 
then  used  to  calculate  theoretical  periods  of  oscillation.  These  theoretical  values 
ulso  differed  in  a uniform  manner  from  the  measured  quantities. 

The  reuson  for  this  difference  is  caught  up  in  the  mechanics  of  the  jig  hangings. 
There  is  a slight  difference  in  the  geometry  when  the  jig  is  hung  from  the  projec- 
tions of  three  vertexes  as  opposed  to  when  it  is  hung  from  projections  of  two 
vertexes  and  an  edge-midpoint. 

Although  the  measurements  performed  on  the  series  of  weights  were  used  to 
calibrate  the  device,  since  only  a single  calibration  constant  was  sought,  the  re- 
sults of  these  measurements  indicate  the  workability  of  the  method, 

The  position  of  the  c.g.  determined  by  the  measurements  was  displaced  about 
I / 2 cm  from  its  true  locution.  Since  the  muss  of  the  subjects  for  which  this  ap- 
proach is  intended  is  on  the  order  of  three  times  greater  than  the  mass  involved 
here,  one  would  expect  that  the  positions  of  the  c.g.  us  calculated  would  be  within 
2 mm  of  their  true  positions. 

The  moment  of  inertia  of  the  series  of  weights  cun  be  visualized  as  a perfect 
sphere.  The  calculated  moment  of  inertia  is  then  a sphere  so  slightly  deformed  that 
its  smallest  and  hugest  diameters  are  each  within  l"„  of  the  true  value. 


Result* 

The  procedures  were  applied  to  nine  subjects,  six  heads,  and  three  head-and- 
necks  at  Tulnne  Medical  School  during  October  5 23,  1971.  The  following  mea- 
surements were  made  for  each  subject: 

1.  Subject  weight. 

2.  Locations  of  right  and  left  auditory  meatuses  in  the  jig. 

3.  Locations  of  right  and  left  orbital  notches  in  the  jig. 

4.  Location  of  anterior  nasal  spine  in  the  jig. 

5.  Weight  distribution  of  projections  for  subject-jig  system. 

6.  Weight  distribution  to  projections  for  empty  jig. 

7.  Periods  of  rotational  oscillation  for  subject-jig  system. 

8.  Periods  of  rotational  oscillation  for  empty  jig. 

These  were  input  to  a computer  program  embodying  the  manipulations  devel- 
oped in  the  mathematics  section  to  solve  for: 

1.  The  c.g.  vector. 

2.  The  moment  of  inertia  tensor. 

3.  The  principal  moments. 

4.  The  alignments  of  these  moments. 

The  value  of  the  measurements  taken  on  the  six  heads  and  the  three  head-and- 
necks  is  questionable.  1 he  weights  of  the  heads  dropped  as  much  as  25“,,  from  the 
time  the  anatomical  sections  were  made  to  the  time  of  the  measurements.  This 
weight  loss  is  due  largely  to  dehydration  in  that  even  though  the  subjects  were  im- 
mersed in  liquid  prescrvat.ve,  the  tissues  tost  their  fluids  and  became  leathery, 

Although  this  dehydration  could  not  be  expected  to  affect  all  tissues  uniformly, 
massless  parameters  of  the  subject  mass  distribution  would  be  less  altered  and 
thus  more  reliable  indicators  of  the  actual  in  vivo  situation  than  those  involving 
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mass.  For  example,  the  c.g.  vector  would  probably  be  more  reliable  than  the 
weight,  The  alignments  of  the  principal  moments  of  inertia  would  probably  be 
more  reliable  than  the  values  of  those  principal  moments. 

For  this  reason,  there  are  tubulated,  along  with  the  principal  moments  of  each 
subject,  the  lengths  of  the  semiuxes  of  an  etlipsoidul  body  of  uniform  density  hav- 
ing the  same  muss  and  moments  of  inertia.  These  parameters  arc  similar  but  not 
identical  to  the  radius  of  gyration  often  calculated  for  bodies  in  rotution  about  a 
single  axis. 

The  data  for  each  subject  include  the  thrce-by*three  moment  of  inertia  tensor, 
the  three  components  of  the  c.g.  vector,  a direction  cosine  matrix  composed  of 
three  row  vectors  of  the  principal  axes  of  the  moment  of  inertia  tensor,  the  three 
principal  moments,  and  the  three  semiaxes  (Fig,  4). 

Although  the  sample  size  in  euch  case  is  too  small  to  vulidute  statistically  any 
strong  conclusions,  the  results  indicute  the  following: 

1 . The  dutu  on  head  3 1 52  are  anomalous. 

2.  The  humun  head  and  heud-and-neck  are  roughiy  symmetrical  about  the  mid- 
sugittal  plane, 

3.  The  anatomical  reference  points  used  to  locate  the  head-referenced  coordi- 
nates are  sufficiently  indicative  of  the  mass  distribution  of  the  head  and  heud-and- 


* - AVERAGED  C.G. 

♦ - ALIGNMENT- DRAWN  THROUGH 

ORIGIN 

t- AVERAGED  ALIGNMENT 
HEAD  NOT  TO  SCALE 


Fig.  5 - Alignments  of  /„  principal  axis  and  c.g.  in  head-referenced  X-Z  (midsaggital) 
plane 
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neck  that  when  these  data  are  described  in  terms  of  these  coordinates,  there  is  a 
noticeable  correlation  between  subjects,  This  is  illustrated  in  Fig,  5,  which  shows 
the  c.g.  and  the  alignments  of  the  i principal  moment  in  the  mldsagittal  plane. 

In  the  future,  lead  cells  wilt  be  substituted  for  the  scule  in  the  c.g.  procedures 
affording  accurate  load  distributions.  A new  hanging  arrangement  will  be  used  in 
the  moment  of  inertia  procedures  allowing  a r .:ch  more  accurate  determination 
of  the  geometries.  These  updated  procedures  will  then  be  applied  to  as  many  sub- 
jects us  possible. 

The  problem  of  weight  losses  In  the  cadaveric  material,  white  probably  not  in- 
surmountable, is  probably  so  dependent  upon  embalming  techniques  that  the  only 
workable  solution  would  be  to  insist  on  the  remains  of  the  recently  departed. 

These  procedures  can  be  readily  adapted  to  other  anatomical  segments.  The 
tetruhedrul  jig  can  be  built  to  the  largest  scale  necessary  without  any  radical 
changes  in  its  structure.  Once  suitable  devices  to  hold  the  subject  segment  within 
the  jig  are  fabricated,  and  the  measurements  and  algorithms  necessary  to  locate 
a three-dimensional  coordinate  system  within  that  segment  formulated,  all  that 
remuins  is  to  apply  the  two  sets  of  procedures  given  here  to  determine  the  c.g,  vec- 
tor and  the  moment  of  inertia  tensor. 


Nomanclatur* 

h„  - nth  component  of  vector  displacement 
C'„,  c„  - nth  constant  in  many  equations  (second  subscript  indicutes  system 
orientation) 

g m acceleration  due  to  gravity  (positive  down) 
t„m  - direction  cosine  between  nth  coordinate  in  one  system  und  mth  coordi- 
nate in  another 
m - system  muss 
p - length  of  pendulum  wires 
T - rotationui  period 
w m weight 

x,  y,  z » cartesian  coordinates 

n.m.f  “ coordinates,  orientations,  or  entries  in  equution 
(..g,  - center  of  gravity 
jig.  j>  J - jig 

.v,  heud  - subject 

sj  - subject-jig  combination 
x - operator  for  vector  cross-product 
T - matrix  transpose 
A » unit  vector 
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Appendix  A 

Thedynumic  equation  of  a balance  scale  (Fig,  A-l)  is 

/ 0 + + (/ 1 +J'i)<i<t>  + (J\ii  - f 1 1, ) - 0 

d « t2 


(A-l) 


where:  •. 

/ - system  inertia 

b m system  dumping 

j\  - weight  being  measured 

J\  « calibrated  weight  used  in  measurement 


Fig.  A-l  - Bal- 
ance scale 


i 


Fig.  A-2  • Pro- 
jection loadings 
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</,  the  offset  of  the  pivot  from  the  lever  arm,  guuruntees  that  the  arm  will  oscil- 
late to  u stable  equilibrium  at  <fi  - 0 whenever,  2 - J\. 

However,  in  the  ease  of  the  c.g.  procedures  (Fig.  A-2),/2  is  not  constant  with  4> 

h - Mh>  + M)/A2 

Oh}  “ — </>  i 2 

d « A 1,^2,  Aj  (A-* 

This  meuns  that  for  small  oscillations,  the  dynamic  equation  of  the  scale  is 

1$  ‘ htj)  ~ (fin  tj/Aj  -/,r/  - fmh\d/h1)<t> 

+ 11/h1  - /|1|)  - 0 (A-: 

Since  the  term  for  in  this  equation  is  negutive,  the  balance  arm  has  no  stable 
position  and  will  ulwuys  move  toward  one  of  its  limits. 

As  an  alternative,  the  test  setup  was  arranged  so  that  the  jig  was  level  when  the 
balance  pan  was  at  its  higher  limit  and  the  balance  arm,  consequently,  at  its 
lower  limit.  The  reading  at  which  this  arrangement  is  no  longer  stable  is  equal 
to  the  load  on  the  bulunoe  pun  times  some  constant  fixed  In  the  scale  geometry. 
This  constant  is 


J‘i  reading  « f 


“ h 


(|t/A  + i^dS 
*2*1  “ <2</A 


where; 

A * angle  of  balance  arm  ut  its  lower  limit 

This  small  error  is  further  alluviated  by  the  least-squares  solution  to  1 2 measure- 
ments distributed  symmetrically  about  the  jig. 

Nomenclature - 

/ - system  Inertia 
b - system  dumping 
d - offset  of  balance  arm  from  its  pivot 
- angular  position  of  balance  arm 
J'iJm  « weight  being  measured 
/ 1 - calibrated  counterweight 
ii,ij  - lengths  of  balance  lever  arms 
I)  = angular  position  of  jig 
/1,,/ij,//.,  - lengths  in  jig  geometry 

• ~ differenthition  with  respect  to  time 
••  - double  differentiation  with  respect  to  time 
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App«ndlx  B 

Given  four  points  in  three-dimensional  space:  (A|,  A),  Z|), . , . , (A4,  A4,  Z4). 
The  midpoint  of  a line  joining  any  pair  of  these  points  is  given  by 


Xij  - (*,  + A,)/2 


similarly 


<»•!) 


where /ty  is  the  midpoint  of  the  line  joining  point  p,  to  point  pj.  Applying  these 

equations,  one  can  determine  pi  j,p2j.Ps4,p«i. 

A further  application  yields  the  midpoints  of  the  lines  joining  pn  to  pi4  and 
p4,  to p2),  and  these  midpoints  are  identical,  being 

(A- 1 + A2  + A3  + A4)/4 

etc.  (B-2) 

etc. 

Thus,  those  two  lines  intersect  so  that  p12,  />2J,  pJ4,  und/>4|  must  becoplanar. 
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Appendix  C 


The  values  lor  t, , and  />,  are  computed  from  the  A',  A,  and  Z positions  of  the 
four  anatomical  landmarks  as  measured  in  the  jig  coordinates. 

Let  (A,,  A,,  Z,)und(A2>  A2,  Z2)  be  the  jig-referenced  positions  of  the  right  and 
left  auditory  meatuses,  respectively,  and  (A.,,  A.,,  Z3)nnd(A4,  A4,  Z4)  are  the  jig- 
referenced  positions  of  the  right  and  left  orbital  notches. 

Then 


and 


b\  - (A)  + A'2)/2 

t>i  - (Ki  + y2)/2 

b>  - (Z,  + Z2)/2 


Mi  - (A’j  + A4  - A,  - A2)/n 

Mi  “ ( Aj  + F4  “ A | - A2)/<» 

Ml  “ (Zj  + Z4  - Z|  - Z2)/« 


where: 


« - KAfj  + a4  - a,  - a2)j  + (Yj  I a4  - a,  - a2)2 

+ (Z4  + Z4  - Z,  - Zi)1)'/'- 

M.i  = t(‘2i(^2  + Z.4  --  kfj  - Z|)  + i,n(A'  + Aj  - A2  - A4))/jd 
Mi  “ 1miIA'2  + A4  - A|  - A.i)  + «,,(/,  + Z|  - Z2  - Z4)]/tf 

MJ  - Imi(A2  + A4  - A,  - A.,)  + i2I(A,  + A,  - A2  - A4)|//i 


(C-I) 


(C-2) 
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(3  - [l.2l(2;  + Z4  - Za  - Z,)  + tji(K|  + Y>  - Y<-  Z4)l2 
+ 1*ji(*s  + *4  - *i  - Afj)  + t|l(Z,  + Za  - Z2  - Z4)]2 

+ u..(y2  + y«  - y.  - yj)  + *2.(*i  + *3  - v2  - *4)i2]'/2 


(C-4) 


■i 

■7)1 


und 


*12  “ *23*31  - *21*33 

‘22  ” *33*u  - *13*31 

*32  " *13*21  - *11*23 


(C-5) 


where: 


*«,.  Y„Z, 


three  components  of  position  oi  point  n in  jig 


b„  - nth  component  of  position  of  origin  of  heud  coordinates 
i nm  ■=  direction  cosine  of  alignment  of  nth  jig  coordinate  to  mth  head 
coordinate 

- quantities  calculated  in  development  for  normalizing  direction 
cosine  but  of  no  permanent  interest 
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Appendix  D 

The  X Y 7.  positions  of  the  four  anatomical  landmarks  in  the  jig  coordinate 
system  are  obtained  with  the  OKI  1760  micromunipulutor(Fig.  D-l),  This  highly 


in  four  ditTcrent  configurations. 

It  is  composed  of  a base  designed  to  attach  to  cither  of  the  two  calibrated  burs 
that  make  up  the  bulk  of  the  DKI  stereotaxic  unit.  The  base  cun  be  moved  along 


Fig.  D-1  • DKI  1760  mlcromanlpu- 
lator 
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where  ti>{  and  <t>2  are  the  ungulur  readings  at  the  first  and  second  pivots. 
In  the  jig  coordinates  for  the  bar  on  the  righi. 


COS  0 I COS  02 


<V 


rc,  + C4  " 

C2  + Ci  - X2 
C’j  + (-b  + X; 


4- 

1 _ 

:>c 

+ 

c c 

1 

C’s  1 

L J 

(O-S) 


For  the  bar  on  the  left, 

0 Oj  |cos<j>|  cos 0} * 
0-10 
0 0 1 


A„- 


c,  + c4 

|C’j  + C 5 - X i 

|£'j  + + ^5 


C'j  + Ad 

(■’ll)  + Q 

r* 


(0-6) 


These  expressions  can  be  simplified  to: 
For  the  bar  on  the  right 


COS  01  COS  02*  ‘ ' 

“(«*  • •" 

c\ 

rc-r  + a? 

. 

• 

• 

• 

(-1  ~ X s 

+ 

c, 

£'\  + X2_ 

(0-7) 


For  the  bar  on  the  left: 

^COS0|  cos  </>;• 


0 
0 

where: 


0 0 

-1  0 

0 1 


'<V 


V, 

C'j  - Aft 

Cl  + x2 


C*4  + Af  I 
C’s  + C'7 


(0-8) 


C„  » lengths  encountered  in  system  hardware 
A',  =■  reading  giving  position  of  manipulator  on  jig 

,Vj  » reading  taken  on  first  arm 
xs  = reading  taken  on  second  arm 

01  = angular  position  about  first  pivot 

02  = angular  position  about  second  pivot 

/(„,  H„  = constants  determined  by  system  configuration 

(’■,  are  constants  determined  by  the  hardware.  A , , 0, , 0>,  A ,,  A , are  too 
readings  for  the  manipulator  in  order  from  the  base  to  the  pointer. 
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Appendix  E 

The  moment  of  inertia  procedures  involve  measurements  of  the  dynamic  be- 
havior of  a trililar  pendulum.  A trifilar  pendulum  consists  of  a rigid  body  sus- 
pended by  three  wires  of  equal  length  attached  to  congruent  points  in  the  body 
and  in  a plane  fixed  perpendicular  to  gravity. 

A complete  statement  of  just  the  kinematics  of  such  a system  would  be  quite 
complex  and  is  beyond  the  scope  of  this  paper.  However,  the  kinematics  for  small 
transverse  and  rotational  displacements  from  the  equilibrium  position  cun  be  ap- 
proximated in  the  following  manner. 

Assume  a rigid  body  suspended  by  three  wires.  Each  of  length  P attached  to 
congruent  points  in  the  body  A,  B,  and  C und  a plunc  fixed  perpendicular  to 
gravity.  The  equilibrium  position  then  has  the  plane  of  A , B , und  C parallel  to  the 
fixed  plane  und  each  of  the  three  wires  parallel  to  the  gravity  vector  if  the  per- 
pendicular from  the  c.g.  to  the  plane  A,  B,C  passes  through  the  triangle  ABC. 

If  the  wires  are  not  permitted  to  go  slack,  this  system  has  three  degrees  of 
freedom.  The  position  and  orientr.ion  of  the  body  in  space  for  small  displace- 
ments (Fig.  E-i)  from  the  equilibrium  condition  cun  be  described  as  functions  of 
the  displacement  i f the  body  in  a plane  perpendicular  to  gravity  and  rotation 
about  un  axis  parallel  to  gravity. 

Assume  a right-handed  coordinate  system  located  in  the  body  with  the  origin 
at  A,  the  T-uxis  passing  through  B , und  the  XY  plane  including  C’such  that  X is 
positive  toward  C. 

Assume  a similar  coordinate  system  located  in  the  laboratory  such  thut  the  lab 
and  body  coordinates  overlap  perfectly  when  the  body  is  in  its  equilibrium 
position. 

Given  small  displacements  of  the  body  relative  to  the  lab  in  x and. a und  small 
rotations  </>  about  the  body  z-axis,  the  displacement  of  A,  B , and  Cin  the  lab 


Fig.  E-l  - Vertical  motion  for 
small  displacements 
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or 


<t>  - — 

6 1 ■*  Cj02  — Ci4>x 

0 2 ~ ‘'2</>2  + 4’2^ 

*c.g.  - * + c4d>  + erf2  + c6<j>y 

yc.t.  - y + erf  + c»<i>2  + c<,<i>x 

z«.g.  - <M*2  + y 2I  + fio</>2  + ‘‘n^  + ‘,i2# 


<v 


where  f|  -C|2  are  all  constants  in  the  system  geometry. 
<-T 


Yti  1 

— Ci  - - 

2 P P 


<■') 


_1_  Yc1  + Xc2  - Y„  Yc 
Ip  Xc 


y«. 


‘’s  “ 2c,gCj 

2a. 
p 

C1  “ YCg 

ZcgYi 


<•* 


2p 


l'V  “ 2/» 


<’io  “ r- 


2/2  L‘* 


Y,.,.  YB  + ( y*  + AT?.  - Fc- 


Ye 


(E-7) 


‘’ii  " “ 


1 c.g. 


<•12 


Y,g. 


(E-8) 


.1 


The  dynamic  equations  for  the  system  can  now  be  developed  from  Hamilton's 
principle  in  the  following  manner.  The  kinetic  coenergy  (TH)  of  the  system  is 

TE  - l/2m|^.g  + A-.g.  + *,2g.) 


+ t/2|/t|<)|  + 2/|2()|02  T 
+ 2/l3<?,  0 + 2/23  0 20  + h)’i>1\ 


and  the  potential  energy  is 


v + mg  zt.  g 


(E-9) 
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forming  the  integral 


(TE  - v)dt 


and  applying  the  variational  calculus  yields 


L 5^. 
[ *■»•  Ox 


, V a2?±  + 2 

+ '**  IT  **  ox 


.1  „ oo , , ..  afl2  , a Ml 

fj  + /,|dl  Tx  + /,li|  IT  + Jl  «T 

..  dfl2  , ,.  M,  , ..  M2 
+ ^21^ i t-  + 'i.^  77  + ^ 77 


02q  O f\ 

+ ms  if  ■ “ 


I U I t)  u^-»-  4.  V 

mI*m  7TT  + ^ “ + 4 


* if]  + ,J'  fy  + /'2*'  ^ + ',2*2  Ty 

' + '»*>  £ + '«*  S’  + '«*  71 


0zc  „ 

+ mg  -7  «*  0 


Jjf  + >>cg  ^ + *.*  + /.,».  f;  + /,2* 

“ /I 


Substituting  from  the  kinemutic  equations, 


f + *•»  if]  + S+  ,i,fi  11+ 

+ /«#  + /J:  + /»<i  ^ + w*  - 0 VB-IO) 


m [Jfc,g,  + Cb<t>j>e,t,  + (2t’y*  + ('ll  0)^c.g,l 

/||‘‘2^l  - h 1^2c2<t> 

- Iuc2<t>4>  + mg[2c9x  + f I i 0 1 - o 
ill  [cb<l>!ic.t.  + ^..g.  + (2c9y  + f|20)*c.g.) 

+ /|2^|t‘j0  + ln^2c2<i>  + ^ 2i*^  f,2 

+ mg\2c9y  + ci2<t>)  - 0 

m [ f c*4  + 2t*j0  + ft.v'V.g,  + (c7  + 2c*0  + 

4-  (2cio0  + t’n*  + ‘•|2y)^.*.l 
4-  "*■  ^12^2  + f|J$  I 1 2<*|  <t>  — n2.xl 

+ 1/22^2  + /|2^1  + /23$](2‘'30  + ^.Vl 
+ /jj<£  + /|3^|  + /2J<?2  + »lgl2(‘l()4>  + ‘'ll-*  + ‘'12^1  “ 0 


(K- 11) 
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Taking  udvantuge  of 


2fVX  + C\\<j>  m i (A-  + C4<j>\  - XCl|, 

2 c9y  + c\i<p  » ^ I.V  + <’7^1  - y».,. 

C 4 Cy 

T|j  « 4*  — C\i  m - — 

1 P P 


+ m | x,.,.  - 0 

+ /«  -p  y*,t,  - o 


C'47«.Vu.,.  + (,4*  + lli<P 

+ m ^,«  - 7 - 7)  * + 7*.*  + ^.*]  -0 


Which  reduce  to 


*«...  + { **.*•  - 0 
J>.*  - { >v„.  ‘ 0 


•^[yya.yjk'.i.y  ll'  y >JL±< 

/„p  r-  * A'c  " *c'r  ' AY 


- X*  - YL  <t>  - 0 


Since  all  the  displacement  terms  are  small,  products  of  these  terms  are  insignifi- 
cant. Eliminating  these  from  the  dynamic  equations  yields 

+ mg[2c,x  + cu<f> ] - 0 
'MJY,.  + mg[2cty  + c,2* ) - 0 
c4mxeii.  + c7My<.4.  + + /ufi,  + /23<>2 

+ wig(2c|0d>  + C||X  + t*i2 y]  - 0 (E-12) 

And  from  the  kinematic  equations 

</>  - — 

0,  & 0 

S;!SO 

xc.,.  « x + erf 

y c.g.  y + cvt> 

z*.,.  **  0 (E- 
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The  equations  in  and  j^g,  describe  the  motion  of  a simple  pendulum  of 
length  p.  The  equation  in  <t>  is  also  a simple  second-order  equation  completely 
uncoupled  from  the  two  transverse  oscillations. 


<t>  (/)  is  a sine  wave  of  period 


- lilt . c 1 ' 
[ m8 


where  c is  u function  of  the  system  geometry 


“ fr'c.g.  Yu 


, x x . x n_y  YhYc  Xi 

■+*  A c.g. A C + A c.g.  v “ A c.g.  v ~ A c.g. 

AC  A C 


This  period  is  the  quuntity  measured  in  the  moment  of  inertia  procedures. 


Nomenclature* 

p - lengths  of  pendulum  wires 
A,  B,  C'  « points  fixed  in  body  geometry 
ifi  « rotation  about  body  2-axis 

0 1 « rotation  about  body  .x-axis 

02  - rotation  about  body  y-axis 

3a  - vector  displacement  in  lab  of  point  a 
x,P  - x and  >>  components  of  that  displacement 
<■„  - constants  in  system  geometry 
TE  - kinetic  coenergy 
m - system  mass 

lHM  » components  of  moment  of  inertiu  tensor 
v m potential  energy 

g m acceleration  due  to  gravity  (positive  down) 

/ - time 

/t,/j  - limits  of  integral 

- single  and  double  differentiation  with  respect  to  time 

^ - partiul  differentiation  with  respect  to  { 

c.g.  - center  of  gravity 
T «*  rotutional  period 
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Appendix  F 


The  computer  program  (Figs,  F- 1 and  F-2)  was  written  for  the  U nivac  1 1 08 
computer  at  the  NASA  facility  in  Slidell,  Louisiana,  No  attempt  was  made  at 
programming  or  formating  elegance, 


(Subject-Jig  Calculations) 
-*-PHIMAS-(c,g,  calculation) 

— R0TIN-(x.  calculation) 

(Empty  Jig  Calculation*) 

— PHIMAS 
— ROTIN 

(Subject  Calculation*) 

—XTRACT- (calculate  tubjtct  c. g,  and  x) 
—HMMNT- (obtain  data  tor  tubjoct  location 
and  alignment) 

•— HDIOR- (colculat*  subject  location  and 
alignmsnt) 

-—TRANS-  (transform  *ubj*ot  parameter* 
to  th*  *ubj*ct  ooordlnat**) 
—VEIGLE- (calculate  principal  moment* 
and  alignments) 


Fig,  F-)  • Main  pro- 
gram 


PHIMAS 

READ  PROJECTION  LOADINGS 
(set  up  matrix) 

[ -CPHI-Uystem  geometry) 

-AXEB- (solve  matrix  equation) 

(back  ealeulat*  projection  loadings) 

[ -CPHI 
‘RETURN 

ROTIN 


rREAD  ROTATIONAL  PERIODS 


(set  up  matrix) 

-RCONS- (equation  constants  from  system 
geometry) 

aeometrv) 


Fig,  F-2  • Phlmas  and 
rotln  programs 


‘•ORNTN-Uyetem  geometry) 

-AXEB 

(bock  calculate  rototlonal  periods) 

-RCONS 

t-ORNTN 
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ADDENDUM 

MAXIMUM  EXPECTED  ERROR 

In  each  of  the  center  of  gravity  measurements,  the  ratio  of  the  measured  projec- 
tion loading,  to  the  total  system  weight,  yields  the  position  of  the  system  c.g.  along 
a distance  Xl,  The  error  In  a single  measurement  AW  would  lead  to  a location  error 
for  the  c.g.  of: 

AX  mM  XL 

CO  WT 

Since  these  c.g.  measurements  are  distributed  symmetrically  about  the  jig,  the 
worst  case  possible  arising  from  comparing  the  twelve  measurements  doubles  the  error. 
Applying  equation  16  leads  to  an  error  In  the  specimen  c.g.  location  of: 

AXcas  a4X  Xt 

This  error  Is  the  magnitude  of  the  worst  possible  vector  dislocation  of  the  subject 
center  of  gravtty  and  assumes  that  the  vector  errors  of  the  two  sets  of  measurements, 
sub|ect-jtg  and  jig,  are  worst  case  and  parallel. 

Errors  In  the  locations  of  the  jig  and  subject-jig  centers  of  gravtty  lead  directly 
to  errors  In  'C',  a function  of  system  geometry  occurring  In  equation  10.  The  worst 
case  error  In  'C'  Is  given  by 

3 AW 

A C 3 2C  R -m-  Xi 


where  R Is  the  distance  of  the  system  c.g.  from  the  geometrical  center  of  the  system 
hanging  points  and  whose  greatest  value  Is  about  5 cm. 

The  resulting  error  to  the  calculated  subject  moment  of  Inertia  Is  about  .5%. 

The  timing  errors  In  the  rotational  period  measurements  are  on  the  order  of  ten 
milliseconds.  Comparing  this  to  a total  measured  time  of  150  seconds  gives  a resulting 
error  In  the  calculated  moment  of  Inertia  of  much  less  than  one  percent. 
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In  aquation  17,  however,  the  calculated  Jig  and  subject-jig  centers  of  gravity  are 
employed  to  calculate  the  subject  moments  of  Inertia: 


From  equation  17 


U. 


i -r  + 1 

U*J  *Uy  g 


r w 
L ,J 


-w 


KM  »J 


KJas 


-W 


KM  J 


Let  W / W « X 

j i 

I -I  - T - 1 x (X  + 1)  W (X  -X  )a 

Us  Hsj  % 9 a kmsj  kmj 

Taking  the  worst  case  vector  errors  in  the  subject-jtg  and  jig  c.g.  locations  as 
parallel  to  the  line  connecting  these  two  points  yields  the  worst  case  error  In  I t 

* 14 

AT  « i 2 A W(2X  + 1)  X,  (X  -X  ) 

M'S  0 ‘•KMSJ  KMJ 

The  second  largest  sourcs  of  error  (the  largest  arising  from  the  deterioration  of  the 
specimens)  then,  Is  that  arising  from  the  errors  In  the  measurements  of  the  projection 
loadings. 

Taking  the  accuracy  of  tl >e  projection  loadtngs  as  .02  pounds,  the  specimen 
wet^t  as  about  8 pounds,  and  the  ratio  of  the  specimen  wei^it  to  the  Jig  weight 
as  3,  yields  as  the  worst  possible  error  for  the  subject  c.g. 


AX 


c»s 


.64  cm 


and  since 


(X  -X  ) < 1,5  cm 

KM  SJ  KMJ 


the  worst  possible  error  for  the  subject  moment  of  Inertia: 

AT,  ■ .40  - 04  KG  -M* 

M'S 

which  Is  less  than  4%  of  an>  of  the  principal  moments  encountered  In  the  specimens. 
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